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1. Introduction 

The low-energy constants (LEC) appearing in the chiral effective Lagrangian, which are of 
great phenomenological importance, can be determined by fitting analytical results from chiral 
random matrix theory (RMT) to lattice data for the eigenvalue spectrum of the Dirac operator. The 
lowest-order LECs are £ and F. While £ can be determined rather easily, e.g., from the distribution 
of the small Dirac eigenvalues, the extraction of F is somewhat more complicated and requires the 
inclusion of a suitable chemical potential [1, 2]. 

Since lattice simulations are restricted to a finite volume, it is important to take into account 
finite-volume corrections to the RMT results, which can be obtained by going to next-to-leading 
order (NLO) in the £ -regime. Recently, finite-volume corrections to the unquenched partition func- 
tion of QCD in the £ -regime were obtained in [3, 4]. However, in order to extract the relevant 
eigenvalue correlation functions the partially quenched partition function of QCD is needed. A 
relatively simple method to obtain the partially quenched theory is to introduce n replicated flavors 
in the unquenched theory and then to analytically continue in the discrete number of quark flavors 
to zero. This so-called replica trick was first used in the theory of disordered systems [5]. It is 
well known that the replica trick is potentially problematic since the analytic continuation from an 
isolated set of points is not uniquely defined. 

In this contribution we choose to use an alternative way to obtain the partially quenched theory 
that does not suffer from the potential problems of the replica trick and can therefore be used 
to check and extend previous results. In addition to the sea quarks, we introduce fermionic and 
bosonic valence quarks. In nuclear physics and condensed matter physics this method is known as 
the supersymmetry method or Efetov method for quenched disorder [6]. In the context of QCD this 
idea was first used by Morel [7]. The effective low-energy theory of QCD with Nf +N V quarks and 
N v bosonic quarks was developed by Bernard and Golterman [8] and by Sharpe and Shoresh [9]. In 
this work we use the effective theory without a singlet particle as discussed by Sharpe and Shoresh 
and consider it in a finite volume and for small quark masses. In order to access F in addition to E, 
we include an imaginary quark chemical potential fL [1,2]. (A first exploratory lattice study of this 
idea was performed in Ref. [10].) We compute the partition function at next-to-leading order in 
the £ -regime and thereby obtain finite-volume corrections of order 1 / \/V to the partially quenched 
theory that translate into finite- volume corrections to the LECs £ and F . Our results agree with 
previous results for the unquenched partition function [3, 4]. The details of this calculation are 
given in a separate publication [11]. 

2. The effective theory 

We define QCD with Nf +N V quarks and N v bosonic quarks by the partition function 

Z = J d[A] <?~ Sym 

where the integral is over all gauge fields A, Sym is the Yang-Mills action, is the Dirac operator, 
m\ , . . . , m^ f are the masses of the sea quarks, m v \ m V N r are the masses of the fermionic valence 
quarks, and m' vl ,. . . ,m' vN are the masses of the bosonic valence quarks. By setting the mass m v ; of 
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a valence quark equal to the mass m' vi of the corresponding bosonic quark, the ratio of determinants 
of this pair cancels and the flavor i is quenched. The corresponding chiral effective theory in the 
partially quenched case of Nf > is defined by the Nambu-Goldstone (NG) manifold given by [9] 



U (x) = exp 



iy/2 



(2.2) 
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where jc(jc) and k(x) are independent N v x (Nf + N v ) matrices with elements in the Grassmann 
algebra, n(x) = 7l(x)^ and 7l'(x) = jv'(xy are traceless Hermitian matrices of dimension N/ + N v 
and N v , respectively, (p(x) G R, and 1„ is the ^-dimensional identity matrix. Note that in the fully 
quenched case of Nf = also the singlet field needs to be included in the effective theory [9]. In this 
contribution, however, we restrict the discussion to the partially quenched case. The Lagrangian of 
the effective theory to leading order in U (x), d p U(x), and the quark mass matrix M is given by 



if = _ Str [d p U(x)- l d p U(x)] - -Str [M*U(x) + U(x)- 1 M] 



(2.4) 



with M = diag(mi, . . . ,m^ f ,m v i, . . . ,m V N r ,m' vl ,. . . ,m[, N ) and low-energy constants £ and F. The 
Lagrangian for nonzero imaginary chemical potential is obtained by replacing the derivative d p 
with the covariant derivative V p defined by 



V p U(x) 



d p U{x) 



iS p0 [C,U(x) 



(2.5) 



where C = diag(jUi , . . . , ji^ f , ji v \ ji v ^ v , ju^ , . . . , ji' vNr ) and ifLj is the imaginary chemical potential 
of quark flavor i. The supertrace Str is a generalization of the ordinary trace that satisfies StrAZ? = 
StrfiA [6]. 



3. Finite-volume corrections to E and F 



In this section we consider the theory in a box of volume V = L0Z4L2L3 in the Euclidean 
formalism. The temporal extent of the box is given by Lq. We use the £-regime power counting 
[12] defined by 



V 



£ . 



(3.1) 



Note that the expansion in e 2 amounts to an expansion in 1 / y/V. To leading order in this power 
counting the fields t, are effectively massless. Therefore, in order to obtain finite propagators, we 
separate the constant mode Uq by the ansatz 



U(x) = Uq exp 



iy/2 



(3.2) 
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with / d 4 x £, (x) = 0. The leading-order Lagrangian is given by 

J^o= l -Str[(d p ^x))(d p ^(x))] -^[lrfUo + U^M] -^-SH[[C,U^)[C,U )) . 
A careful analysis of the propagator, see Ref. [11], yields 



8ad8bc{-l) eb ~ TT^abdcd 
N f 



where 



fd[Z}e-J d4 *^o ' 
A(x) is the massless propagator without zero modes [13], and 
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fori <b<N f +N v , 
^ 1 for Nf + N v <b< N f + 2N V . 



(3.6) 



Note that the propagator does not depend on the valence quark number N v . The propagator A(x) is 
finite in dimensional regularization and contains the dependence on the volume. The contributions 
to the Lagrangian at next-to-leading order in e are given by 



Sf 2 = ^2+^2+^2 



with 
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We will integrate out the fluctuations in E, in order to obtain an effective finite-volume partition 
function. The term Jzf 2 M couples to Uq and M and corrects the leading-order mass term to 
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where A(0) in dimensional regularization is given by [13] 



A(0) 



(3.12) 



The so-called shape coefficient j3i only depends on the quantities /,• = Lj/V 1 / 4 with i = 0, 1,2,3. 
The terms in J*f 2 c couple to Uq and C and correct the leading-order chemical potential term to 
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with 



d 4 x (doA(x)Y 



2VV 



&00 



where &oo is another shape coefficient [14] that only depends on Iq, . 
effective low-energy constants £ e ff and F e g, 



Nj-l 
NfF 2 



F 2 



A(0), 



A(0)- / d 4 x (d A(x)y 



(3.14) 

, I3. Thus we can read off 

(3.15) 
(3.16) 



This is the same result as previously derived for the unquenched theory [3, 4]. The terms in 5£J£ 
and also potential contributions from the integration measure are not relevant for finite-volume 
corrections to £ and F at next-to-leading order in e, see Ref. [11] for details. 



4. The constant-mode integral 

In the infinite- volume limit the fluctuations in t, are suppressed and the theory becomes zero- 
dimensional. It is therefore described by RMT, and the integral over the constant mode Uq defined 
in Eq. (3.2) should recover the RMT result. A proof for the case of vanishing chemical potential 
is given in Ref. [11]. We also note that fixing topology to a single topological sector v merely 
amounts to a change of the integration manifold of Uo and therefore does not alter the discussion 
of the finite-volume corrections given above. For a detailed discussion of the theory in a single 
topological sector we refer to Ref. [11]. 



5. Conclusions 



In this work we have calculated the partially quenched partition function of QCD at next-to- 
leading order in the e-expansion at nonzero imaginary chemical potential. We considered a theory 
with Nf+N v fermionic quarks and N v bosonic quarks, as formulated by Sharpe and Shoresh [9], 
in a finite volume V with microscopic quark masses M, i.e., MVL = 0(e). The knowledge of the 
analytic form of the partially quenched partition function suffices to obtain all spectral correlation 
functions of the Dirac operator lj>. In this sense our results for the finite-volume behavior of the 
theory hold universally for all observables that can be obtained from spectral correlation functions 
of ift. We found that the partially quenched partition function has the same finite- volume corrections 
as the unquenched partition function of QCD with Nf quarks. 

In Fig. 1 we show the finite-volume corrections at NLO to the low-energy constants £ and 
F as a function of the box size L in a symmetric box. Note that the effects of the finite volume 
increase with the number of sea quark flavors Nf and that, depending on Nf, a box size of 3 — 5 
fm is necessary to reduce the effects of the finite volume at NLO to about 10%. The effects are 
calculated at F = 90 MeV. In Fig. 2 we show the effect of an asymmetric box with Nf = 2 and 
L = 2 fm. An important message of this figure is that the magnitude of the finite- volume corrections 
can be significantly reduced by choosing one large spatial dimension instead of a large temporal 
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Figure 2: Effect of an asymmetric box with parameters Nf = 2, L = 2 fm, and F = 90 MeV. We compare 
a large temporal dimension Lq with L\ = Li = L3 = L (left) to a large spatial dimension L3 with Lq = L\ = 
L 2 = L (right). 



dimension. The reason for this behavior is that the chemical potential only affects the temporal 
direction, see Eq. (2.5), and therefore breaks the permutation symmetry of the four dimensions. 
This manifests itself in the propagator 

f d A x (d A{x)) 2 (5.1) 

which, as shown in Eq. (3.14), contains a term proportional to Lq/\/V, where Lq is the size of 
the temporal dimension. This term leads to an enhancement of the corrections in case of a large 
temporal dimension. Choosing instead one large spatial dimension, the finite-volume corrections 
are reduced, unless the asymmetry is too large. For the parameters used in Fig. 2, the optimal value 
is L3/L w 2. 

This is good news. Many lattice simulations (at zero chemical potential) are performed with 
L\ = L2 = L3 = L and Lo = 2L. To determine F, it suffices to introduce the imaginary chemical 
potential in the valence sector. Therefore, one can take a suitable set of existing dynamical con- 
figurations and redefine Lq <-> L3 before adding the chemical potential. 1 This will minimize the 
finite-volume corrections for both £ and F, at least for the parameter values chosen in Fig. 2. 

'Note that this procedure increases the temperature of the system by a factor of two. One needs to check that the 
system does not end up in the chirally restored phase, in which our results no longer apply. 
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